a_angle = (
            (
                self.m # NOTE: THIS WAS PREVIOUS SQUARED WHEN IT SHOULDNT BE
                * self.einstein_radius ** (self.slope - 1.0)
                * radial_grid ** (2.0 - self.slope)
            )
            / (self.m**2.0 - (3.0 - self.slope))
            * self.k_m
            * np.sin(self.m * (polar_angle_grid - self.angle_m))
        )